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^^ ABSTRACT 

Z^ Discerning the radiative dissipation mechanism for prompt emission in Gamma- Ray Bursts (GRBs) 

requires detailed spectroscopic modeling that straddles the uFy peak in the 100 keV - 1 MeV range. 

r^ Historically, empirical fits such as the popular Band function have been employed with considerable 

success in interpreting the observations. While extrapolations of the Band parameters can provide 

some physical insight into the emission mechanisms responsible for GRBs, these inferences do not 

provide a unique way of discerning between models. By fitting physical models directly this degeneracy 

can be broken, eliminating the need for empirical functions; our analysis here offers a first step in this 

direction. One of the oldest, and leading, theoretical ideas for the production of the prompt signal is 

PQ the synchrotron shock model (SSM). Here we explore the applicability of this model to a bright Fermi 

HH GBM burst with a simple temporal structure, GRB 090820A. Our investigation implements, for the 

^; first time, thermal and non-thermal synchrotron emissivities in the RMFIT forward- folding spectral 

r^ analysis software often used in GBM burst studies. We find that these synchrotron emissivities, 

0^ together with a blackbody shape, provide at least as good a match with the data as the Band GRB 

spectral fitting function. This success is achieved in both time-integrated and time-resolved spectral 

fits. 

Subject headings: acceleration of particles — gamma-ray bursts: individual (GRB 090820A) — gamma 
^ rays: stars — methods: data analysis — radiation mechanisms: non-thermal — 

radiation mechanisms: thermal 
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1. INTRODUCTION 

In the most popular paradigm for gamma-ray bursts 
of both long and short durations, it is typically assumed 
that prompt 7-ray emission results from the dissipation 
of kinetic energy in a relativistically expanding fireball 
mediated by multiple internal shocks (e.g. see Piran 
1999, or Meszaros 2001, for reviews). These shocks are 
presumed to diffusively accelerate a fraction of the elec- 
trons from thermal upstream distributions to higher en- 
ergies. Usually only particles in the exponential tail of 
the Maxwellian are available for acceleration. Thus, for 
relativistic shocks, the expected outcome is that the par- 
ticle distribution consists of a Maxwellian with a power- 
law tail at high energies. Based on this scenario, the 
radiative emission should consist of two components, 
quasi-thermal and non-thermal photons from electrons 
spiralling along magnetic field lines in optically thin re- 
gions of the jet. There could also be an additional pho- 
tosphcric contribution of Planckian form, originating in 
distinct, optically thick environs, perhaps interior to the 
regions spawning synchrotron emission. 

To date, the characterization of GRB spectra has been 
dominated by the use of the empirical Band function 
(Band et al.JIT993), a parametrized, smoothly broken 



power law that was devised in the era of the BATSE 
experiment on the Compton Gamma-Ray Observatory 
(CGRO). Several authors have used measurements of the 
Band spectral shape parameters to infer properties of the 
physics involved in GRB emission. In particular, the fit- 
ted spectral indices defined by N^ (x E^" below the vF^ 
peak and iV^ ex E^p above it, may be compared with val- 
ues predicted from synchrotron emission: the low-energy 
self-absorption index, a, of -1-1 (in photon flux units), 
the synchrotron 'line of death' index of —2/3, the 'sec- 
ond line of death' at the fast cooling value of —3/2, the 
high energy index, /3, characterizing power-law particle 
acceleration, and the various spectral differences between 



these ([Preece et al.|1998l|Lloyd fc Petrosianl200Cil [Lloyd 
Ronning & Petrosian|2U02l|Preece et al.|2U02|). However 



Ronning fe Petrosian|2u62| " Preece et al.|2U02[ ). However, 
it becomes difficult to discern between models through 
the Band function when the fitted low energy indices rep- 
resent a power law only asymptotically, and when many 
models predict similar Band indices. In fact, the Band 
function's inherent shape and curvature only loosely ap- 
proximates the shape of the applicable physical models 
making it difficult to draw conclusions about emission 
mechanisms directly from Band function fits. A way to 
break this degeneracy is by fitting more realistic emission 
models to the data, which in addition provides deeper in- 
sights into the burst environment. In Section [2] we detail 
the emission model that we use to fit GRB spectra. We 
present this model as a first step. In future work we will 
explore additional models in an attempt to discern be- 
tween them. We describe our observational results with 
this model in Section [3l 

2. MODEL AND MOTIVATION 

We propose to test an emission model composed of 
synchrotron emission and a thermal blackbody. This 
model is the most general form of the standard fireball 



model. Non-thermal synchrotron emission is historically 
the most favored process invoked to explain prompt GRB 
signals. The motivation for the inclusio n of a blackbody 
component comes not only from theory ( Goodman [1986 



Meszaros3'2002'; Roes & Meszaros '2005) and prevfous 



searches ( Ryde and Pe'cr ,,2009) but also the recent dis- 
covery o f a significant quasi-thermal compone nt in GRB 



spectra ( [Guiriec et al.|2010 Ryde et al.|2010 ). However, 
in Guiriec et al. (20lfi; GRB 1007^4liJ as well as Ryde 
et al. (2010; GRB 090902B), the non-thermal portion of 
the spectra is approximated by empirical functions that 
lack direct associations with the physical parameters. 

In order to model optically-thin synchrotron emission 
in a physical way, we adopt the parametrization pre- 
sented in Baring & Braby (2004, hereafter BB04), which 
was modified slightly from the choice of Tavani (1996). 
Theory and numerical simulations predict that the elec- 
tron energy distribution resulting from diffusive shock ac- 
celeration should be composed of two components (e.g. 
see Baring 2011 for an overview), which to first order 
can be approximated by a superposition of a relativistic 
Maxwellian and a super-thermal power-law tail; 



rieh) = "0 



J-7/7T J^ 



ii) 
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where Q{x) is a step function with Q(x) = 1 for a; > 1 
and zero otherwise, and 7t is a measure of the post- 
shock electron temperature. This is a quasi-isotropic 
distribution, in the comoving frame of reference of the 
GRB outfiow (the mildly-relativistic speed of an internal 
shock in this frame does not change this form signifi- 
cantly), with the dependence on pitch angle being omit- 
ted for simplicity, though it can be incorporated in the 
no factor. The shock acceleration electron distribution 
therefore depends on five parameters, three of which, the 
power-law index, S , the relative normalization, e (which 
can be related to the acceleration efficiency), and the 
product ry7T , which defines the minimum Lorentz factor 
of the power-law, pass unmodified into the expression of 
the photon flux and are thus flt parameters for the GRB 
data. In Tavani's original exposition 77 was fixed to unity 
and the power law component smoothly joined to the ex- 
ponential portion of the Maxwellian (i.e., with virtually 
no discontinuity). This would be the case of 'saturated' 
acceleration, where all of the electrons above the peak in 
the Maxwellian have been accelerated. BB04 indicated 
that values 77 ~ 3 and e ;$ 0.1 closely reflect populations 
usually found in simulations of shock acceleration even 
ones based on diverse and contrasting approaches (e.g. 
see Niemiec & Ostrowski 2004; Spitkovsky 2008; Baring 
2011; and references therein). For simplicity and general 
facility of spectral fitting, we adopt the compact form in 
Eq. (n]), deferring direct fitting with specific simulation 
modeToutput to future studies. Here, 77 = 3 is adopted 
as a representative value that incurs no significant dis- 
continuity in transitioning from the Maxwellian to the 
non-thermal population when e ;$ 0.1 . 

In a truly physical model, the electron distribution 
function should be perfectly continuous, contrasting 
Eq. (II]). Here, we have left both rj and e as parame- 
ters free to vary, observing that folding the distribution 
with the synchrotron emissivity function in Eq. ([2]) below 
yields continuous emission spectra. Thus, while not ex- 
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phcitly joining the two components of Eq. (II]) smoothly, 
the subsequent fitting of GRB spectral data provides a 
robust and informative indication concerning the relative 
contribution of each component, as was done in BB04. 
More precise modeling with truly continuous electron dis- 
tributions is left for future investigations, but is unlikely 
to alter the essential conclusions of our work here. 

To determine the radiation flux, F^{£) ex £N^ , emit- 
ted by these electrons, this distribution is convolved with 
the standard synchrotron emissivity (e.g. Rybicki and 
Lightman, 1979; see also BB04): 



/OO 



where 



^ ( f ) d7 



J^ [w) — w K^/^{x)dx 



(2) 



(3) 



expresses the single-particle synchrotron emissivity (i.e. 
energy per unit time per unit volume) in dimensionless 
functional form. The characteristic scale for the syn- 
chrotron photon energy is 



c i B „ 2 2 



(4) 



where B^-c — 4.41 x 10 Gauss is the quantum critical 
field. When convolved with the distribution in Eq. (fTl), 
the substitution 7 — )■ tj^t- in Eq. Q then defines the 
scale for the break energy of the synchrotron continuum 
resulting from the truncated power-law portion of the 
distribution (see Table 1). In modeling prompt burst 
emission, the relativistic nature of the outflow introduces 
an extra parameter, the bulk Lorentz factor F of the 
flow, which blueshifts the spectrum so as to introduce 
the r factor in Eq. Q, so that Eq. (l2| then expresses 
the synchrotron flux m the observer's frame. Accord- 
ingly, while the electron distribution parameters 5 and 
e can be constrained by prompt emission spectroscopy, 
the precise values of 7t and the environmental quanti- 
ties B and F are indeterminate, being subsumed in the 
single parameter Trj^^^B that is defined by a spectral 
fit in a given time interval. 

For fits where non-thermal synchrotron components 
dominate, the energy of the vFi, peak determines the 
value of the peak energy. Well below this structure the 
flux index is -1-1/3 and well above it, the flux index is 
— ((5 — l)/2 . This is the simplest synchrotron model to 
consider. Strong cooling synchrotron models possess a 
similar mathematical character, but elicit a gentler break 
and a steeper spectrum below the break that is often 
more difficult to fit to observations. Treatment of such 
cooling models, and inverse Compton scenarios will be 
deferred to future work. We note also that models where 
e ^ 1 and the non-thermal synchrotron component is 
small or insignificant, the high energy tail of the ther- 
mal synchrotron component is necessarily exponentially 
declining with energy. Such forms have severe difficulty 
in fitting GRB spectra that possess extended power-law 
tails, a common occurrence, yielding e ^ 0.1 as an antic- 
ipated frequent inference in this GRB spectroscopy pro- 
tocol. 

In summation, our emission model consists of a two- 
component synchrotron function (thermal and power- 



law), plus a blackbody, all boosted from the outfiow 
frame, by the bulk Lorentz factor F , to the observer's 
frame. Along with the blackbody component, this spec- 
tral model has seven fit parameters; values for two of 
these parameters, 77 and 5 , are fixed for reasons de- 
tailed in Section [3] Owing to the intensive numerical 
integration involved, such functions have previously not 
been used for forward-folding spectral fitting, particu- 
larly in the CGRO/BATSE era. We have implemented 
this photon model into the RMFIT spectral analysis soft- 
ware and demonstrate our technique by fitting fitting the 
prompt emission of GRB 090820A; one of the brightest 
GBM bursts with simple temporal structure. 

3. OBSERVATIONS 

On 20 August 2009, at To=00:38:16.19 UT, the 
Gamma-Ray Burst Monitor (GBM) onboard the Fermi 
Gamma-ray Sp ace Telescope trigge red on the very bright 
GRB 090820A ( |Connaughton|2009|). This GRB also trig - 
gered Coronas Fhoton-Rl'-2 ('C hakrabarti et al ||2009[ ). 
The burst location was initially not in the F'OV of the 
Large Area Telescope (LAT) onboard Fermi but was 
bright enough to result in a Fermi spacecraft repointing 
maneuver. However, Earth avoidance constraints pre- 
vented such a maneuver until 3100 sec after the burst 
trigger and the burst was not detected at higher energies 
by the LAT. The most precise position for the direction 
of the burst comes from the GBM trigger data which lo- 
calizes the burst to a patch of sky centered on RA = 87.7 
degree and Dec = 27.0 degree (J2000) with a 4 degree 
error, statistical and systematic. The current best model 
for systematic errors is 2.8 deg rees with 70% weigh t and 
8.4 degrees with 30% weight ( jBriggs et ar]|2011[ ). We 
verified that our analysis does not change significantly 
using instrument response functions for assumed source 
locations throughout this region of uncertainty. 

GBM is composed of 12 sodium iodide (Nal) detectors 
covering an energy range from 8 keV to 1 MeV and two 
bismuth germanate (B GO) detectors sensiti ve between 
200 keV and 40 MeV (|Meegan et al.||2009|). Figure |l] 
(top two panels) shows the light curve of GRB 09082uK 
as seen by GBM, from 8 to 200 keV in the Nal detectors 
(top) and from 200 keV to 40 MeV in the BGO detector 
(bottom). GBM triggered on a weak precursor which we 
do not include in the analysis. The main light curve be- 
gins at To -I- 28.1s. The main structure of the light curve 
consists of a fast rising pulse with an exponential decay 
lasting until To-l-60 s. A second, less intense, peak begin- 
ning at To-|-30 s is superimposed on the main peak. With 
such a high intensity and simple structure, this GRB 
allows for detailed time-resolved spectroscopy. Because 
this burst is intense, calibration issues make the Iodine 
K-edge (33 keV) prominent in the count spectra owing 
to small statistical uncertainties, and we remove energy 
channels contributing to this feature from our spectral 
fits. In addition, an effective area correction is applied 
between each of the Nal detectors and the BGO during 
the fit process. This correction of « 23% is used to ac- 
count for possible imperfections in the response models 
of the two detector types. 

We simultaneously fit the spectral data of the Nal de- 
tectors with a source angle less than 60 degrees (Nal 1 
and 5) and the data from the brightest BGO detector 
(BGO 0) using the analysis package RMFIT. We use a 
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forward- folding technique that convolves the detectors' 
response with the proposed photon model to generate a 
count spectrum to compare to the data; the parameters 
of the photon model are then adjusted so as to optimize 
the Castor C-stat statistic. The Castor C-stat differs 
from Poisson likelihood by an offset which is a constant 
for a particular dataset. 

We perform a fit to the integrated spectrum and 
find that it is best represented by synchrotron emission 
from thermal and power-law distributed electrons with 
an additional blackbody component characterized by a 
kT w 42 keV (C-Stat/DOF = 558/353). The vF^ spec- 
trum is displayed in Figure [2] and the best-fit values in 
Table [T] We also performed a fit using the Band function 
(C-Stat/DOF = 593/355). We find in concordance with 
BB04 that emission from power-law synchrotron dwarfs 
the emission from thermal synchrotron by at least 3 or- 
ders of magnitude. The value of 77 is fixed to 3, the choice 
adopted by BB04: it is a value that accommodates distri- 
butions typically determined by shock acceleration simu- 
lations. When fitting the power-law synchrotron compo- 
nent we have to fix the value of the power-law index to its 
best fit value to remove a correlation between the ampli- 
tude and the index; this does not change the fit statistic 
but does mean that the amplitudes obtained are valid 
only for that index. The inferred electron distribution 
from this fit is shown in Figure |4] We note that the in- 
ability to simultaneously constrain the power-law index 
and amplitude of the synchrotron function may be solved 
in future studies by including joint fits with LAT data, 
whenever available. 

For the time resolved analysis we fit four bins labeled 
a, b, c and d as shown in Figure |T] with the various 
synchrotron models. The corresponding electron distri- 
butions inferred from these fits are displayed in Figure 
tWe also fit the Band function to each spectrum to 
ow that in nearly all cases the physical models can fit 
the data as well as the Band function. We chose the time 
binning by finding a balance between high signal-to-noise 
and evolution of the spectral shape so that we can iden- 
tify the time evolution of each component throughout the 
burst. Where possible, we fit all three components simul- 
taneously. Due to the similarity in the spectral shapes 
of the low energy portions of the thermal synchrotron 
and power-law synchrotron components it is not always 
possible to constrain all of the fit parameters especially 
when one component is much stronger than the other. 
Therefore, when one component is dominant we include 
only that component in the fit. The ability to fit both 
components in the time integrated fit is most likely due 
to the fact that both components are significant over the 
interval. 

From bins b to c the spectrum is best described by syn- 
chrotron emission from power-law distributed electrons 
in addition to a blackbody (Table IT] and Figure l4|. The 
thermal synchrotron component is too weak to mean- 
ingfully include it in the fit. We find that the inten- 
sity of the power-law synchrotron increases significantly 
from bin b to c while the blackbody component remains 
nearly constant in intensity. The spectral index of the 
electrons in these intervals varies from -4.4 to -5.9. Such 
values are consistent with those expected from diffusive 
acceleration theory, for the specific case of superluminal 



shocks ( [Baring 2011[ ), i.e. those where the mean mag- 
netic field angle to the shock normal is significant. This 
geometrical requirement establishes efficient convection 
of particles downstream of relativistic shocks, thereby 
steepening their acceleration distribution. The black- 
body component decreases in intensity at this point but 
the temperature remains constant within errors. In bins 
a and d, with weaker emission, several models are es- 
sentially statistically tied. It is possible that PLS-(-BB 
persists throughout the entire GRB. Alternatively, the 
GRB could even begin in bin a with thermal synchrotron 
emission and transition to the PLS+BB emission. If this 
were true we would be seeing emission from electrons 
that have not yet been accelerated into a power-law dis- 
tribution by the shock. The C-stat values for all of the 
models fit in each bin are displayed in Table [2] 

While it is not possible to constrain all parameters in 
all the bins, it should be stressed that this is due to 
natural correlations in the synchrotron functions. These 
difficulties do not arise when using the Band function 
because it has a simpler parametrization. 

4. DISCUSSION AND CONCLUSIONS 

In this paper, we have shown that thermal and non- 
thermal synchrotron photon models, with an additional 
blackbody, are well consistent with the emission spec- 
tra of GRB 090820A in various time intervals. These 
are physical models that afford the ability to constrain 
parameters that are physically meaningful, for example 
key descriptors of the electron distribution that is mo- 
tivated by shock acceleration theory. By implementing 
these models into a forward-folding spectral analysis soft- 
ware we have been able to directly constrain many of the 
physical model parameters and their respective errors; a 
first in the field of GRB spectroscopy. This constitutes 
substantial progress over the use of the empirical Band 
function to fit prompt GRB spectra, which has been a 
nearly universal practice to date. The results presented 
here enable more rigorous statements about the validity 
of GRB emission models, moving the study of prompt 
burst emission into a new era. 

Our modeling has focused on the standard synchrotron 
shock model with the addition of a blackbody compo- 
nent. The spectral fitting reveals a complex temporal 
evolution of the separate components. While spectral 
evolution is a well-known feature of GRBs, this type of 
fitting can enable direct physical interpretation of the 
evolution. These fits provide evidence that the line of 
death issue (Preece et al. 1998, 2002) can be overcome 
naturally with a combination of synchrotron and black- 
body emission: the prominence of a blackbody com- 
ponent with its flat Rayleigh- Jeans portion would de- 
rive a comparably-fitted Band function with a fiat low- 
energy index. This was also suggested by Guiriec et al. 
(2010) where the authors used simultaneous fits of the 
Band function and a blackbody. Note that it is possible 
that other physical models may, in fact, produce supe- 
rior fits to the data for GRB 090820k and other bursts. 
Strongly-cooled synchrotron emission, inverse Compton 
and jitter radiation are popular candidates, and our work 
here motivates the future development of RMFIT soft- 
ware modules for these processes. 

A principal finding of the analysis in this paper is that 
the power-law synchrotron component is orders of magni- 
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tude more intense than the thermal synchrotron compo- 
nent during the peak of the burst, the latter contribut- 
ing at most a few percent of the flux. This confirms 
the finding of BB04 for BATSE/EGRET bursts GRB 
910503, GRB 910601 and GRB 910814, which was a 
theoretically-based perspective that did not fold mod- 
els through the detector response matrices. They had 
noted that full plasma and Monte Carlo diffusion sim- 
ulations of shock acceleration clearly predict a power- 
law tail in the particle distribution that smoothly ex- 
tends from the dominant thermal population (e.g. see 
also Baring 2011, and references therein). This tail is 
several orders of magnitude smaller than what is found 
when fitting synchrotron emission to burst spectra. It 
is not clear how such non-thermally-dominated distribu- 
tions can arise near shocks, providing a conundrum for 
the standard synchrotron shock model. Limited smooth- 
ing of the sharp peak of the non-thermal electron com- 
ponent will not alter this conclusion. 

This result is also in accord with Guiriec et al. (2010), 
in their analysis of GRB 100724Q, who fitted its GBM 
spectra with a combination of the Band model and a 
blackbody. They too found that an unrealistically high 
efficiency for the acceleration mechanism or a source size 
smaller than the innermost stable orbit of a black hole 
was required to invoke the standard fireball model for 
explaining the origin of the 7-ray emission. Therefore, it 
was surmised therein that the outflow from the jet was 
at least partially magnetized. 

To conclude, the success of this analysis in isolating 
the relative contributions of a handful of distinct spec- 
tral components indicates that it is imperative for the 
field of GRB spectroscopy to move away from the use 
of the empirical fitting functions: many physical mod- 
els can asymptotically approximate the Band spectral 
indices, rendering it difficult to discern between them 
particularly near the vFi, peak. Instead, direct compar- 
isons of the fitted physical models are possible, and are 
required to truly discriminate between the various emis- 
sion processes. The fitting of physical SSM/blackbody 
spectra here offers a clear advance beyond empirical fits. 



and provides the impetus for further development and 
deployment of physical modeling of prompt burst emis- 
sion spectra. 
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Figure 1. Light curve of GRB 090820 A as observed by GBM. Ttie two panels show the count rate in the two NAI detectors (top) and 
BGO (bottom). The dashed lines indicate the time intervals (a, b, c, d) used for the time-resolved analysis (see Figure 3 and Table 1). It 
is clear that the burst consists of two main peaks and that this burst is very bright in the BGO detectors. 



Table 1 

The fit parameters for the time-integrated (first row) and time-resolved spectra. The fit parameters for the blackbody component are its 

amplitude (Agg) and energy (kT). The fit parameters for the non-thermal components are described in section[2] The break energy 

£t = £c(7 ~> 'IIt) corresponds to employing the substitution 7 — > tj^t in Eq. (J4|l. Note that the ratio of the amplitudes is not equal to 

the ratio of the fluxes. 



Time interval 



Model 



"-0 (7'' ^cm ^keV ^) 



£c (keV) 



i^s^'^cm'^keV''^) kT (keV) 



Time integrated 


TS-I-PLS+BB 


3437+0-^1)4 


8711-1 


10 3o+U-^M 
iU.jy_g 245 


a 




TS 


Q7O+0.189 
^•■^'^-0.176 




».dt)i_Qg3 


b 




PLS-(-BB 


859t«st.? 


— 


14 9/I+0-848 


c 




PLS-fBB 


I.901IHI3 X 10* 


- 


15.22lO«^ 


d 




TS-I-BB 


2.196tS:If6 


- 


^■035t°fil 



4.9=" 3.0 



4,4a 



3.0 
3.0 



1.774 
1.818 



-1-0.410 
-0.356 
-1-0.400 



X 10 

X 10 



-0.344 

8.383+t?? X 10" 



42.27 



-1-1.49 



35.32ti?? 



QQ '7+2.13 

3»-'_1.92 
28.40+?Iq 



fixed at best-fit value 



Constraints on the Synchrotron Shock Model 



Times: 28.672: 48.129 s 



1000 



> 



100 



00 



5 



-5 
-10 



BGO_00 
NALOl 
NAL05 




+ 



**H|^^^ 4l*#^^rt„^^ 



10 



100 



1000 
Energy (keV) 



10000 



Figure 2. The integrated spectrum of GRB 090820A. We are able to resolve three components, thermal synchrotron, power-law syn- 
chrotron, and a blackbody. Energy channels near the Nal K-edge are omitted. The deviations in the fit residuals are the due to systematics 
in the detector response resulting from the high count rate and spectral hardness of this burst. However, deviations are never greater than 
4(T and do not significantly impact the values of the best fit parameters. The multiple curves near the peak of the spectrum are an artifact 
of the effective-area correction applied to each detector and not related to the different fitted models. 



Time Interval 



Band TS 

464/355 466/357 
432/355 742/357 



TS + BE PLS 
464/355 467/357 
445/355 555/357 



450/355 1088/357 488/355 558/357 
404/355 421/357 403/355 406/357 



PLS + BE 
465/355 
434/355 
434/355 
405/355 



Table 2 

The c-stat per degree of freedom for each time model in the selected time intervals. 
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Figure 3. The electron distribution corresponding to the integrated spectrum. The non-physical jump in the amplitude between the 
Maxwellian and the power-law distribution (parametrized by e) at r) is clearly seen. 
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Figure 4. The time-resolved spectra for GRB 090820A. The spectra represent bin a with thermal synchrotron only (top left panel), bin 
b with power-law synchrotron -|- blackbody (top right panel), bin c again with power-law synchrotron + blackbody (bottom left panel), 
and finally bin d with thermal synchrotron -|- blackbody (bottom right panel). As with Fig. p] the multiple curves are associated with the 
effective area correction. 
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Figure 5. The electron distributions for the tiinc-rcsolvcd spectra. The choice of rj with a power-law only distribution is arbitrary due to 
the fact that £c and rj both scale Ep^^ak- 



